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Abstract 

We present a pilot-wave model for quantum field theory in which the Dirac sea 
is taken seriously. The model ascribes particle trajectories to all the fermions, in- 
cluding the fermions filling the Dirac sea. The model is deterministic and applies 
to the regime in which fermion number is superselected. This work is a further 
elaboration of work by Colin, in which a Dirac sea pilot-wave model is presented 
for quantum electrodynamics. We extend his work to non-electromagnetic inter- 
actions, we discuss a cut-off regularization of the pilot-wave model and study how 
it reproduces the standard quantum predictions. The Dirac sea pilot-wave model 
can be seen as a possible continuum generalization of a lattice model by Bell. It 
can also be seen as a development and generalization of the ideas by Bohm, Hiley 
and Kaloyerou, who also suggested the use of the Dirac sea for the development 
of a pilot- wave model for quantum electrodynamics. 

1 Introduction 

In the pilot-wave interpretation of deBroglie and Bohm for non-relativistic quantum 
theory [1-3], also known as Bohmian mechanics, systems are described by both their 
wavefunction and by particle positions. The particles move along deterministic trajec- 
tories, under the influence of the wavefunction. The incorporation of particles in the 
description of quantum systems makes it possible to provide a description of reality in 
which rather ambiguous notions such as 'measurement' and 'observer' play no funda- 
mental role. 

In this paper we present a way of formulating a pilot-wave theory for quantum fields. 
In order to put the model into context, we first briefly review the history and status of 
pilot-wave theories for quantum field theory. 

The history of extending non-relativistic pilot-wave theory to quantum fields starts 
with Bohm's seminal paper [3]. In that paper Bohm gives a pilot-wave interpretation 
for the quantized electromagnetic field. In this interpretation the additional variables, 
also called beables, which together with the quantum state make up the description 
of the quantum field, are not particle positions. Instead the quantum field is also 
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described by an actual field configuration. In a similar way one can construct a pilot- 
wave interpretation for various other bosonic quantum fields, like for example the scalar 
field, see [4-17]. 

In contrast to what might be expected, the introduction of field beables for fermionic 
quantum fields does not appear to be so straightforward. There is an attempt by 
Holland [12,18], who took a field of angular variables as the beable and an attempt 
by Valentini [7,8], who took anti-commuting fields as beables. However, in the case 
of Holland's model some questions need to be addressed before it can be considered 
a viable model [17]. Valentini's pilot- wave model on the other hand is much more 
problematic. Although it might seem natural to introduce anti-commuting fields as 
beables for fermionic quantum fields, it is unclear how this could be done in a way such 
that the standard quantum predictions are reproduced [16, 17]. 

Inspired by these difficulties, Struyve and Westman [19] proposed not to introduce 
beables at all for fermionic quantum fields, only beables for the electromagnetic field. 
This forms a minimalist pilot-wave interpretation, which is nevertheless capable of re- 
producing the standard quantum predictions. 

A different approach consists of introducing particle positions as beables. This ap- 
proach has appeared to be more successful for fermionic quantum fields. Bell presented 
a model for quantum field theory on a lattice [20], where the beables are the fermion 
numbers at each lattice point. Bell's model differs from the usual pilot-wave program 
in the fact that it is indeterministic. However, Bell expected that the indeterminism 
would disappear in the continuum limit. 

Sudbery [21] was the first to consider a continuum generalization of Bell's model. 
He considered a non-relativistic particle, without the possibility of particle creation or 
annihilation, and found that the equation of motion for the particle reduced to the one 
originally proposed by de Broglie and Bohm, which is deterministic. This result was later 
also obtained by Vink [22] . Much later, two different possible continuum generalizations 
were presented that were capable of dealing with particle creation or annihilation. There 
is the work by Diirr et al. [23-27] in which a stochastic generalization is presented and 
the work by Colin [28-30] which deals with a deterministic model. 1 

The fact that there are two generalizations of Bell's lattice model for the continuum 
originates in a different reading of Bell's work. Diirr et al. take the fermion number to be 
the number of particles plus the number of anti-particles. Their model therefore ascribes 
trajectories to both particles and anti-particles. The trajectories are deterministic until 
there is particle creation or annihilation. 

Colin has a different reading of fermion number. In standard quantum field theory 
the 'fermion number' is often taken to mean the number of particles minus the number 
of anti-particles. We call this number F. Instead of using the particle-anti-particle 
picture one can also use the Dirac sea picture, in which an anti-particle is a hole in 
the sea of negative energy particles. In the Dirac sea picture of quantum field theory, 
which is generally considered to be equivalent to the particle-anti-particle picture, there 

1 Bell's original model presupposes that space is discrete whereas time is continuous. The possibility 
of generalizing Bell's model for discrete time has been studied in [31]. 
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is a natural notion of the fermion number, namely the number of positive plus negative 
energy particles. We call this number F^. This number just equals F plus the number 
of particles in the Dirac sea, the latter being an infinite constant. We keep calling both 
F and Fd the fermion number. Colin now takes the fermion number to mean Fd. 

So, in the work by Colin, which deals with quantum electrodynamics, the picture of a 
Dirac sea is taken seriously. Position beables are introduced for the particles in the Dirac 
sea, i.e. the negative energy particles, and for the positive energy particles. It is the fact 
that the fermion number Fd, as the number of positive plus negative energy particles, 
is always positive, together with the fact that it can be considered superselected in 
quantum electrodynamics, that allows the construction of a deterministic model in a 
similar fashion as in non-relativistic pilot-wave theory. Fermion number superselection 
means that the quantum state cannot evolve into a superposition of states with different 
fermion number. In particular superselection of fermion number implies conservation of 
fermion number, i.e. that the fermion number operator commutes with the Hamiltonian. 
(Note that fermion number conservation does not imply fermion number superselection. 
Although conservation of fermion number implies that the fermion number of a quantum 
state will never change during Schrodinger evolution, it could still change during collapse, 
in which case the fermion number would not be superselected.) 

The model that Colin found as a continuum version of Bell's model for lattice quan- 
tum field theory is a model that was suggested earlier, though not developed in detail, 
by Bohm, Hiley and Kaloyerou, see [10] and [11, p. 276]. In this paper, we further 
elaborate on this model. We study to what extent the model can be generalized as to 
account for other types of interactions, such as the weak and strong interaction. In 
addition, we study how the pilot-wave model can be regularized, i.e. how it can be made 
mathematically well defined. The regularization is required at two levels, namely at the 
level of the Schrodinger equation which describes the evolution of the quantum state, 
and at the level of the equations describing the beables. Instead of introducing a lattice 
regularization, like Bell, we introduce a cut-off regularization. This allows us to keep 
continuous space and to preserve the determinism of the pilot-wave model. We also 
study how the pilot-wave model reproduces the predictions of standard quantum field 
theory. 

As mentioned before, a key ingredient in the construction of the Dirac sea pilot-wave 
model of Colin and Bohm et al. is that the fermion number can be considered supers- 
elected in the case of electromagnetic interactions. In particular, the fermion number 
is conserved. To date, no experiment has shown a violation of fermion number con- 
servation, not even if non-electromagnetic interactions are involved [32]. Moreover, all 
Feynman diagrammatic processes in the Standard Model preserve the fermion number. 
A familiar example is the weak process d — > u + e + i> e , which takes place in /3-decay. 
Nevertheless the Standard Model does predict a violation of fermion number conserva- 
tion. This predicted violation of the fermion number is a non-perturbative effect, i.e. 
it does not appear in perturbative treatments, such as Feynman diagrammatic ones. 
The violation is a prediction of the electroweak theory [33] and only appears for suffi- 
ciently high energies. According to one recent estimate [34], a violation is expected to 
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be observed no sooner than in the Very Large Hadron Collider (the hypothetical future 
hadron collider with performance significantly beyond the Large Hadron Collider). 

Our generalization of the Dirac sea pilot-wave model of Colin and Bohm et al. will 
hold in the regime where the fermion number is conserved. Although the pilot-wave 
model can therefore not reproduce all the predictions of the Standard Model, it can at 
least reproduce the predictions of the Standard Model in some low energy regime. For 
example the pilot-wave model can perfectly describe the regime in which the electro- 
weak symmetry is broken by spontaneous symmetry breaking. In this paper, we will 
therefore restrict our attention to this energy regime. The model of Diirr et al. does 
not rely on the conservation of any quantum number and therefore seems capable of 
describing fermion number violating interactions. 

Next to fermion number superselection, also the introduction of a Dirac sea is a 
key ingredient in the construction of the pilot-wave model. In a filled Dirac sea all 
negative energy states, corresponding to quarks and leptons with all possible quantum 
numbers, are occupied. Superselection and positivity of the fermion number Fd then 
imply that we can introduce position beables in a similar way as in non-relativistic 
quantum mechanics. The position beables are further not distinguished by any label 
corresponding to properties such as mass, spin, charge or flavour. This is similar to 
Bell's model for lattice quantum field theory, since the beables in this model are the 
numbers of fermions at each lattice point. 

With the introduction of position beables also for the negative energy fermions, the 
Dirac sea pilot-wave model yields an ontology which is not immediately related to our 
image of a classical world. Furthermore, in standard quantum field theory measurements 
are not expressed in terms of fermion number densities which include the negative energy 
fermion contributions. As a result, it is not so straightforward to show that the pilot- 
wave model reproduces the predictions of standard quantum field theory. Instead of 
providing conclusive evidence, we only give arguments that support our expectation 
that the model reproduces the standard quantum predictions. 

We also want to mention that the pilot-wave model is not Lorentz covariant at the 
fundamental level. However, insofar as the model reproduces the predictions of standard 
quantum field theory, Lorentz covariance is regained at the statistical level. There exist 
some attempts to construct a Lorentz covariant pilot-wave model, see [35] for a good 
review. On the other hand, see also [36] for an argument that the fundamental symmetry 
of pilot-wave theories is Aristotelian rather than Lorentzian. 

The outline of the paper is as follows. First, in Section 2 we describe in general the 
technique by which the pilot-wave model is constructed. Then in Section 3, the actual 
model is presented in the context of the effective quantum field theory which arises in 
the Standard Model after spontaneous symmetry breaking. In Section 4, we discuss a 
cut-off regularization and indicate how the model reproduces the standard quantum field 
theoretical predictions. Finally, in Section 5 we add a note on the predicted violation 
of fermion number in the standard model. 
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2 General framework 



Before we pass to the general framework for constructing the pilot-wave model, we 
first recall the pilot-wave interpretation for non-relativistic quantum theory that was 
originally proposed by deBroglie and Bohm [1-3]. In this interpretation, a quantum 
system is described by its wavefunction ^(x 1; . . . , x n ,t) in configuration space M. 3n for 
all times t, which satisfies the non-relativistic Schrodinger equation 

<9?/;(xi,...,x n ,t) ( ^h 2 V\ \ 
ih — = I - 2^ 7^- + V{X lt . . . , X n ) I ^(Xi, . . . , X n , t) , (1) 

and by n particles, which move in physical space M 3 , and for which the possible trajec- 
tories Xfc(t) are solutions to the guidance equations 

" V * = ^T' (2) 

where ip = expiS/h. 

How the pilot-wave theory reproduces the standard predictions of quantum theory 
has been discussed before in many places, see e.g. [3]. A key ingredient is that the 
positions of the particles are distributed according to \ip\ 2 over ensembles where all the 
systems have wavefunction ip. The special distribution \ip\ 2 is called the equilibrium 
distribution [37,38]. The particular form of the guidance equations guarantees that if 
the distribution of particle positions equals IV'Cxlj • • • > x n, ^o)| 2 at a certain time to, then 
it will equal IV^Xi) • • • ; x n ,t)| 2 at times t. This means that the distribution keeps its 
functional form of ip. This property is also called equivariance [38]. Possible justifications 
of quantum equilibrium have been discussed in [7, 11, 37, 38]. 

There is a simple way to derive the pilot-wave interpretation. Our Dirac sea pilot- 
wave model for quantum field theory will be derived along similar lines. This derivation 
goes as follows. In standard non-relativistic quantum theory |^>(xx, . . . , x^, t)\ 2 is in- 
terpreted as the probability density to find the particles in the respective volumes d 3 Xi 
around the positions x,, in a measurement on an ensemble of systems all described by 
the wavefunction ip. The conservation of probability yields the continuity equation 

f + EVr(»Jltf)=0, (3) 

fc=l 

which can be derived from the Schrodinger equation (1). If one insists that particles 
have positions at all times, i.e. also if no measurement is performed, then the analogy 
with fluid dynamics suggests the equation of motion (2) for the particles. In this way 
we arrive at the pilot-wave interpretation. 

Note that originally, neither deBroglie nor Bohm followed the above derivation to 
arrive at the pilot- wave theory [39] . De Broglie postulated a precursor of the guidance 
equation before the wave equation was known, so that the guidance equation was not 
obtained from considering the wave equation first. In his Ph.D. thesis in 1924 deBroglie 
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namely postulated that the momenta of particles should be proportional to the gradient 
of a phase, where the phase corresponds to an unspecified periodic phenomenon. Only 
later, in 1926, Schrodinger provided a wave equation for this phase. On the other hand, 
when Bohm rediscovered the pilot- wave theory in 1951, he was in the first place mo- 
tivated by the resemblance to classical Hamilton- Jacobi theory, which arises when the 
Schrodinger equation is rewritten in terms of the polar decomposition of the wavefunc- 
tion. 

Nevertheless, this way of deriving the pilot-wave theory provides a scheme which can 
be used (and is often used) to construct pilot-wave models for other quantum theories. 
Let us sketch this scheme. We restrict ourselves to the case in which particle positions 
are introduced as beables. The scheme is immediately amenable for the introduction of 
other type of beables, like fields, strings, etc. 

Suppose we have a Schrodinger equation 



ih 



d\ m) 

dt 



H\m) 



(4) 



with Hamiltonian H and a positive operator valued measure (POVM) P(d?x\ . . . d 3 x n ) 
for position. This means that (ip(t)\ P(d 3 x\ . . . d 3 x n ) \ip(t)) gives the probability to find 
n particles in the respective volumes d 3 Xi around the positions Xj in a measurement 
at time t. 2 This probability determines a probability density p^(xi, . . . ,x n ,t) by the 
relation 

p*(x 1 ,...,x n ,t)d 3 x 1 ...d 3 x n = (^(t)\P(d 3 x 1 ...d 3 x n ) \iP(t)) . (6) 

The Schrodinger equation implies the following time evolution for the probability den- 
sity: 



dp^xi, . . . ,x n ,t) 
dt 



d 3 xi...d 3 x n + (ip{t)\ih P(d 3 Xl ...d 3 x n ),H \if>(t)) = . (7) 



If we can write the second term in the form 



n 

(ij(t) | Ik [P(d 3 X 1 . . . d 3 X n ), H\ [V(*)> = Vfc ' ( V M X 1> • • - X n, *)) d 3 Xl . . . d 3 X n , 

k=l 

where the are some vector fields, then we can construct a pilot-wave model with 
particle positions as beables, by postulating the guidance equations 



'k ■ 



(9) 



The equilibrium density for the particle positions is given by . 



2 For example, in the case of non-relativistic quantum theory we have a position PVM, i.e. a projection 
valued measure, which is a special case of a POVM. The position PVM is given by 

P(d 3 x 1 . . . d 3 x n ) — |xi, . . . ,x„) (xi, . . . , x„| d?xi . . . d 3 x n , 

where the states |xi, . . . ,x„) are the standard position eigenstates. 



(5) 
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Note that we can not always write the second term in (7) in the form (8), for example 
when the number of particles is not conserved. In such cases one still has the pilot- wave- 
type models of Diirr et al. which include stochastic jumps [23-27]. 

3 The pilot-wave model 

In this section, we present our pilot-wave model for quantum field theory. We start with 
the Hamiltonian of the Standard Model after spontaneous symmetry breaking, so that 
the weak interaction bosons and the fermions have acquired their mass. This means that 
we do not have to worry about the predicted violation of fermion number conservation, 
which only appears for higher energies. We further discuss the predicted violation of 
fermion number in Section 5. 

At this stage, we do not introduce a regularization of the theory. We postpone the 
introduction of a regularization to Section 4. 

In this and the following sections we use units in which h — c—1. 

3.1 The Hamiltonian 

The Hamiltonian is of the form 

H = H B + H F + Hi . (10) 

The part H B is the part of the Hamiltonian that contains only bosonic fields, i.e. the 
interaction fields and the Higgs field. It contains the free Hamiltonian of the bosonic 
fields together with all the coupling terms between the various bosonic fields. The part 
Hq is the free Hamiltonian for the fermionic fields, i.e. the quark and lepton fields. Hj 
contains the interaction terms between the fermionic and bosonic fields. 

The Hamiltonian depends on the Dirac field operators Va(x) of the quarks and 
leptons. The label A comprises quantum numbers like lepton number, quark number, 
quark colour, etc. 3 The Dirac fields ip and their conjugates ^ satisfy the equal-time 
anti-commutation relations 

{ft,„(x), fi W (x')} = 5 xy 8 aa ,8(x - x') . (11) 

The other fundamental anti-commutation relations are zero. The labels a and a' are 
Dirac spinor indices. We give a representation for these operators in the following 
subsection. 

For the bosonic operators we do not introduce an explicit quantization scheme. We 
leave the corresponding bosonic Hilbert space and its orthonormal basis {|£)} unspec- 
ified. We expect that the subtleties in quantizing the bosonic fields do not affect the 

3 Note that we take the neutrinos to be Dirac particles. The possibility still exists that neutrinos 
are Majorana particles. In that case the neutrino would be its own anti-particle, which would lead to 
a violation of fermion number conservation. 
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key properties of our pilot-wave model, at least not in the presently considered energy 
regime. 

The free Hamiltonian for the fermionic fields Hq and the interaction Hamiltonian 
Hi are respectively given by 

Ho = S d 3 x$l(x)[-ia ■ V + /3m A ]^ A (x) , (12) 

A "* 

Hl= Yl /^L( X )^W,a'(x)^V,a'(x), (13) 
A,A',o,a' 

where m A is the mass of the particle with quantum numbers A. The kernel /i A ,a,A',a'( x ) 
carries a hat since it involves the bosonic field operators. It contains no fermionic 
field operators. The form of the interaction Hamiltonian is typical for an interaction 
introduced by minimally coupling the free fermionic Hamiltonian Hq with gauge fields. 
Also the Yukawa interaction which couples the fermionic fields to the Higgs boson is of 
this form. 

Instead of writing out the bosonic Hamiltonian H B or the interaction Hamiltonian Hi 
explicitly, let us just present parts of the latter here. For the electromagnetic interaction 
Hamiltonian, for example, the kernel is given by 

(x) = qx5 xy (7VL, A„(x) , (14) 

where q\ is the charge of the particle with quantum numbers A and the vector 
potential. Another example is the interaction Hamiltonian which represents the coupling 
of the leptons to the weak interaction vector bosons Wr^. In order to have a simple 
presentation of this Hamiltonian let A stand for just the lepton flavour. We write 
A = (l,s), where I represents the lepton family, i.e. electron, muon or tau family, and 
where s = 1,2 indicates the electrically charged or neutral lepton (i.e. the neutrino) in 
the family. With this notation the kernel reads 

hx*x'A*) = l^pfiv (7V(1 - 7 5 )) Qa , (Wyi%(x) + Wrt^(x)) , (15) 

where g is the coupling constant for the group 577(2) in the electroweak theory. 



3.2 The Dirac sea picture 

The standard Fourier expansion of the Dirac field operator is given by 



?a(x) = £ 



d 3 p 



c A , s (p)n A , s (p)e 4p - x + dl>K s (p)e-^ x 



(16) 



where -E A (p) = \/p ■ P + m \- The operators c i Xs and c\ s respectively create and an- 
nihilate particles with quantum numbers A and helicity s. The operators cft x s and d\, s 
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respectively create and annihilate the corresponding anti-particles. They satisfy the 
relations 



{c AiS (p),4' lS '(p')} = hx'S ss 'S(p--p') , {£,«(p),dA',«'(p')} = h\'S S s'S(p-p') ■ (17) 

The remaining fundamental anti-commutation relations are zero. 

Instead of expressing the Dirac field operators in terms of creation and annihilation 
operators for particles and anti-particles, which have positive energy, we can express 
them in terms of creation and annihilation operators for particles of positive and negative 
energy 

feW = E / ^§yf y^pT K*(pK*(p) eip ' x + Ca,.(pK.(-pK px ] • (is) 

The operators (\ s and (\ >s respectively create and annihilate particles with negative 
energy. They are connected to the anti-particle operators by the relations C\,s(p) = 

2U-p)- 

In standard quantum field theory the state containing no particles or anti-particles 
is noted as |0). In the absence of any interactions this state would be the vacuum. In 
the Dirac sea picture this state is obtained from the Dirac sea vacuum \0d), which is 
the state that contains no particles of positive or negative energy, by applying all the 
creation operators of negative energy particles, i.e. 

i°> = n sl.(p)i°^> ■ ( i9 ) 

A,s,p 

In this picture, an anti-particle corresponds to a hole in the Dirac sea. 

Historically, the Dirac sea picture was convenient to interpret the Dirac equation 
as a many particle equation that was able to explain the appearance of pair creation 
in quantum electrodynamics. For non-electromagnetic interaction the Dirac sea picture 
was abandoned in favour of standard quantum field theory [40]. In standard quantum 
field theory particles and anti-particles can be freely created, without the need to refer 
to the Dirac sea picture. However, in order to construct our pilot-wave model it will 
appear advantageous to reintroduce the Dirac sea picture. It will allow us to formulate 
a pilot-wave theory in which particles travel along deterministic trajectories. 

3.3 The fermion number operator and position PVM 

The fermion number operator is given by 



F d = J d 6 xF d (x) , (20) 

with 

F d (x) = £^ a (x)fe, a (x) (21) 

A, a 
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the fermion number density operator. The fact that fermion number indeed corresponds 
to the total number of positive and negative energy particles, or equivalently the total 
number of particles minus the number of anti-particles plus an infinite constant, can be 
seen from the expression 




= W d 3 p(^ s (p)c AiS (p) + d AiS (p)4 )S (p)) 

X,s •* 

= E / ^ (4>Fa iS (p) - + W) • ( 22 ) 

X,s 

The state with zero fermion number is the state \0d)- 
Let us now introduce the states 

|xi,Ai,ai;...;x n ,A n ,a n ) = (x x ) . . . ^ A (x n )|0 D ) . (23) 

\Jn\ 

These states are normalized as 

( x i) Ai, a\] . . .] x n , A n , a n |x' 1; A' 1; a[; . . . ; x' ra , A^, a n ) 

= ~7 Z^ -1 )^*^ 1 ) ~ *'l) S K W K 6 ^iA ■ ■ ■ 5 ( X -W - 0$X a{n) X'Ja a{n) a> n , (24) 

(7 

where the sum runs over the permutations of the indices 1, . . . ,n and where (— 1) CT is 
the sign of the permutation. Making use of the commutation relations 

= $, a (y)<5(x - y) , (25) 

one can easily show that the states |xi, Ai, a\\ . . . ; x„, X n , a n ) are eigenstates of the 
fermion number density operator 

n 

^d(x) |xi, Ai, a x \ ■ ■ ■ ; x n , A n , a n ) = ^ <5(x - x 4 ) |x x , A x , a x ; . . . ; x n , A n , a n ) . (26) 

i=l 

The corresponding eigenvalues are degenerate since they do not depend on the Aj or aj. 

It is clear that the states |x 1; Ai, a\\ . . . ; x n , A n , a n ) represent, for all values of A, and 
Oj, states for which there is a fermion at the positions xi to x n . The total fermion 
number is equal to n. Therefore the PVM 

P(d 3 £i . . . d 3 x n ) = ^2 |xi,Ai,ai;...;x n ,A n ,a n ) ® |0 

x <g> (xi, Ai, ai; . . . ; x n , A n , a n | rf 3 xi . . . d 3 x n (27) 



F d (y),g ia (x) 



10 



can be interpreted as a position PVM (remember that the states |£) form an orthonor- 
mal basis of the bosonic Hilbert space). The fact that the measure P(d 3 X\ . . . d 3 x n ) is 
projection operator valued is straightforwardly checked by using the normalization (24) 
of the states |xi, Ai, ai, . . . ; x n , A n , a n ). 

The position PVM P{d 3 x\ . . . d 3 x n ) also corresponds to the conventional position 
operator in the Dirac theory when the Dirac sea picture is adopted, see e.g. [41]. It is for 
example for this position operator that Schrodinger derived his famous Zitterbewegung 
[42]. Note that the position eigenstates |xi, Ai, ax; . . . ; x n , A n , a n ) are superpositions 
of positive and negative energy states. Therefore these states are sometimes called 
'abstract' position eigenstates, to contrast them with states which are built out only of 
positive energy states, like for example the eigenstates of the Newton- Wigner position 
operator [41]. 

Now we would like to regard the states |xi, Ai, a\\ . . . ; x n , A n , a n ) with n G N as a 
basis for our Hilbert space. However, the Hilbert space spanned by these states would 
not contain states of physical interest. Physically relevant states, which are states which 
contain a finite number of positive energy particles and a finite number of holes in the 
negative energy sea, namely have an infinite fermion number. However, dealing with 
an infinite number of particles would lead to a number of complications in maintain- 
ing mathematical rigour. First of all we would have the field theoretical divergences. 
Secondly, since there does not exist a generalization of the Lebesgue measure to an in- 
finite dimensional space, one cannot assume the naive infinite particle number limit in 
expressions like those for the position PVM [35]. We will come back to these issues in 
Section 3.6. For the moment, we keep working with states with finite fermion number. 
The hope is that some key properties of the pilot-wave model survive when passing to 
some regularized theory. In Section 4, we discuss in detail such a regularized theory, 
which is obtained by introducing an ultra-violet momentum cut-off and by assuming a 
finite space. 

3.4 Fermion number superselection 

Since the fermion number operator commutes with the Hamiltonian H, the fermion 
number is conserved. We further assume that fermion number is superselected. As 
explained in the introduction this assumption is empirically adequate, although that 
may change in the future if the predictions of Standard Model turn out to be correct. 
The fermion number superselection implies that we can restrict the fermionic Hilbert 
space, which is spanned by the states |x 1; Ai, a±] . . . ; x n , X n , a n ) with n e N, to a subspace 
spanned by states |xi, Ai, a±] . . . ; x n , X n , a n ), where the fermion number n is fixed. 



In the Hilbert space with fixed fermion number n, we can expand a quantum state 
\jj;(t)) as 




x. 



n 



i,ai,...,A Tl ,a„,£V x l' • • • 



x n, t) |xi, Ai, ai] . . . ; 



x n , A n ,a„) ® |f) , (28) 
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where the expansion coefficients ipx 1) a 1 ,...,x n ,a n ,$( yi ij • • • >* n ,t), which are given by 

V'Ai,oi,...,An,on,f( x i> • • • , x„, t) = ((x l7 Ai, ai; . . . ; x n , A n , a n | <g> (f |) , (29) 

are anti-symmetric under permutations of the indices 1, . . . ,n. In terms of these ex- 
pansion coefficients, the position density that corresponds to the position PVM (27) is 
given by, cf. equation (6), 

/(xi,...,X n ,t) = \^\ 1 ,a 1 ,...,\ n ,a n ,i(xi,---,Xn,t)\ 2 . (30) 

Ai,ai,...,A n ,a n ,^ 

Since we have a fixed fermion number, the fermion number density operator can be 
written as 

/n 
J25(x-x l )P(d 3 x 1 ...d 3 x n ). (31) 
i=l 

For later purposes, we also introduce the operator 

F d (B) = [ d 3 xF d (x) , (32) 

JB 

which is the operator corresponding to the number of fermions in the region B. The 
associated expectation value (ip(t)\Fd{B) \ip{t)) gives the average number of fermions 
in the region B. This expectation value can also be expressed in terms of the position 
density p^ as follows: 

m)\F d {B)\m)= I d s x [y25(x-x l )(m\p(d 3 x 1 ...d 3 x n )\m) 

Jb J i=1 

= d 3 x d 3 xi... \ d 3 x n ^ £(x - x^)p v '(xi, . . . , x n , t) 
Jb Jri Jr3 



= n d 3 Xl / d 3 x 2 ... d 6 x n pV(xx, . . . , x„, t) , (33) 
Jb Jr 3 Jr 3 

where in the last equality we have used the fact that the density is completely 
symmetric in its arguments. 

3.5 Pilot-wave model 

We are now ready to present our pilot-wave model. As mentioned before we assume a 
finite fermion number. We will further elaborate on this issue in Section 3.6. 

We introduce position beables for all the fermions, including the fermions in the Dirac 
sea. For the equilibrium density we take the density p^ which corresponds to the position 
PVM P{d 3 X\ . . . d 3 x n ). As explained in Section 2, we can find guidance equations by 
considering the continuity equation for this density. First note that it follows from 
equation (7) that the parts of the Hamiltonian that commute with the position POVM 
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do not contribute to the guidance equation. Since in our case the bosonic Hamiltonian 
H B and interaction Hamiltonian Hj commute with the position PVM P(d 3 X\ . . . d 3 x n ), 
we have that only Hq contributes to the guidance equations. The fact that the bosonic 
Hamiltonian H B commutes with the position PVM is straightforward. The proof that 
the interaction Hamiltonian Hj commutes with the position PVM is given in Appendix 
A. 

With the only contributions coming from the free fermionic Hamiltonian Hq, we 
find that the velocity field is given by 

V t= Yl ^A 1> o 1 ,..,A n ,a n ^a£l... > o n; o i ,... l <V'A 1 ,ai,...,A„ l < > f/p^, (34) 

Al,...,A„,£ 

<»!■■••> »n 

where 

a a 1 ,...,a n ;a' 1 ,...,a' n = ^ a i a 'i ' ' ' a a k a' k ■ ■ ■ &a n a' n ■ (35) 

The guidance equations for the particles are now defined to be 



dt 



< ■ (36) 



These guidance equations are the natural generalizations of the guidance equations for 
the many particle Dirac equation [11]. The dynamics given by (34)-(36) is also a special 
case of a dynamics proposed in [23-27]: it is the 'minimal process' associated with the 
position PVM P(d 3 xi . . . d 3 x n ) and the Hamiltonian H, see e.g. Section 5 in [27]. 

The equilibrium distribution is given by . Having an equilibrium distribution is 
a key requirement for showing that the pilot-wave model reproduces the predictions of 
standard quantum theory. We discuss this in more detail in Section 4.3 in the context 
of a regularized model. 

Because the wavefunction V'Ai,ai,...,A n ,a„,£( x i> • • • , x n,^) is anti-symmetric under per- 
mutations of the label 1, . . . , n of the particle positions, also the particle dynamics, given 
by (34)-(36), is invariant under such permutations. Therefore the particles are identical, 
i.e. they are not distinguished by any property such as mass, charge, flavour which could 
be associated to the label. 

Since the particle label is actually redundant (it has no physical meaning), one could 
formulate a pilot-wave dynamics on the configuration space M? n /S n , where S n is the 
group of permutations a of the set {1, . . . ,n} [43], instead of on the configurations space 
M. 3n . Alternatively, one could formulate a pilot-wave dynamics on the configuration 
space that is obtained from the space M. 3n /S n by removal of configurations for which 
two or more particle positions are the same [43]. This latter configuration space can 
also be identified with the space n M. 3 = {S C IR 3 !^^ = n}, which consists of n-element 
subsets of M. 3 [44]. The dynamics on M. 3n /S n or n M 3 can be obtained from the dynamics 
on R 3n by a suitable projection. 4 



4 Interestingly, a dynamics on "R 3 can actually be formulated in general, i.e. also for quantum 
theories in which systems are not necessarily described by completely anti-symmetric wavcfunctions 
[44]. 
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In Bell's pilot-wave type model on the lattice, the beables are the fermion numbers 
at each lattice point [20]. Hence, this approach also makes manifest that particles are 
considered to be identical. In our model, we can introduce a fermion number density 

n 

F d (x) = ^5(x- Xi ) (37) 
i=i 

corresponding to the fermion positions xi, . . . ,x n . The fermion number density is the 
natural generalization of Bell's fermion numbers to the continuum. In the rest of the 
paper, we interchangeably use the positions or the number density to characterize the 
fermions. 5 

The fermion number density should not be confused with the fermionic (electric) 
charge density, 6 for which the operator is given by 

g(x) = ^g A ^ a (x)^ AiQ (x), (38) 

X,a 

where q\ is the charge of the particle with quantum numbers A. The charge operator Q = 
J d 3 xQ(x) represents the sum of the charges of the positive and negative energy particles. 
The states |xi, Ai, a\\ . . . ; x n , A n , a n ) are also eigenstates of the fermionic charge density 
operator, with eigenvalues Y^i=\ ^^(x — X j). 

When Colin presented the Dirac sea pilot-wave model for quantum electrodynamics 
[28-30], only electrons and positrons were considered, in which case the fermion number 
density operator is proportional to the fermionic charge density operator. So in this 
particular case, one could equally well say that beables are introduced for charge density. 
In the model presented in this paper this can not be said; the beables correspond solely 
to fermion number density. 

Fermion number superselection played a key role in the construction of the pilot- wave 
model. The pilot-wave model as such obtained now implies an even stronger notion of 
fermion number superselection, which is not available in standard quantum theory, and 
which is called strong superselection. This is explained in detail in [45]. So note that if 
there was no superselection according to standard quantum theory, then our pilot-wave 
model would fail to reproduce the standard quantum predictions, simply because our 
model implies this even stronger notion of superselection. 

3.6 The need for regularization 

As mentioned before, the theory needs regularization. The regularization is required 
because of the usual field theoretical divergences. Further, in order to have well defined 

5 Note that we attached several distinct meanings to the notion of particle. We have the quantum 
mechanical notion of a particle, whereby a particle is represented by its quantum state, and the notion 
of a particle in pilot-wave theory, where it is described by its position. The quantum mechanical notion 
of a particle can further be understood as in the standard representation of quantum field theory, which 
uses 'particles' and 'anti-particles', or as in the Dirac see representation. A similar ambiguity arises in 
our use of the notion of 'fermion'. Wc hope it is clear from the context which notion we have in mind. 

6 The total charge density consists of the fermionic charge density plus the bosonic charge density. 
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probability densities for the position beables, we need to deal with the non-existence of 
a generalization of the Lebesgue measure to an infinite dimensional space. 

A possible way to regularize the pilot-wave model is as follows. We could assume 
that space is bounded and impose some suitable boundary conditions. For example we 
could assume that space is a cubic box of volume V and that we have periodic boundary 
conditions for the field operators. For the Fourier expansions of the field operators ^a( x ) 
this means that we have to replace 



d 3 p 




HE' (39) 



where the possible momenta in the sum on the right hand side are given by p = 
(ni,n2,n 3 )2ir/ L, where (711,72,2,713) G 1? and L is the length of a side of the cubic 
box. This makes the number of possible momenta countable. By further assuming an 
ultra-violet momentum cut-off we have that the number of possible momenta is finite. 
This implies that any possible state can only contain a finite number of particles. In 
particular, the filled Dirac sea will only contain a finite number of particles. There- 
fore also in the pilot-wave model we only need to introduce a finite number of position 
beables. We discuss this way of regularization in more detail in the following section. 
The introduction of this regularization will cause a slight modification of the equations 
presented in the previous section. 

Apart from regularization problems, there would be another problem when we would 
take the naive limit of the number of fermions to infinity in the pilot-wave model pre- 
sented in the previous subsection. The problem is the identification of an image of a 
familiar macroscopic world in the beable configuration [35]. We would expect that we 
just have to count the number of fermions in different regions in physical space in order 
to find such an image. However, this procedure would not work for the following reason. 
We have seen that the quantity (ip{t) \ Fd(B) \i/j(t)) gives the average number of fermions 
in the region B. For states of physical interest, i.e. states which describe a finite number 
of positive energy particles and a finite number of holes in the negative energy sea, the 
average number of fermions in every region of non-zero measure, is infinite. For example, 
for the state |0) we have (see Appendix B for the details of the calculation) 

(0\F d (B)\0) = ^ 3 J2 J d 3 P J >* = oo. (40) 

This means that every region in physical space looks the same if only the number of 
fermions is counted in those regions. As such, an image of the macroscopic world would 
clearly be absent. It is an open question whether there exists an alternative way of 
identifying such an image. The problem disappears when an ultra-violet momentum 
cut-off is introduced, since then 

ci Um i»> = ^ £ / p|sA *pj B ** = ^ ? /, *' < 00 ' (41) 

We will discuss this in more detail in the following section. 
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4 Ultra-violet momentum cut-off regularization 



In this section, we regularize the theory by imposing an ultra-violet momentum cut-off 
|p| < A. As indicated in Section 3.6, imposing the cut-off is actually not sufficient to 
completely regularize the pilot-wave model. We mentioned that a complete regulariza- 
tion could be obtained by assuming a bounded space and by assuming some suitable 
boundary conditions. We will not take this last step in regularizing the theory, since it 
would merely complicate some of the calculations, without changing the qualitative con- 
tent. Nevertheless, at some places in this section we will implicitly assume a bounded 
space (e.g. in Section 4.3 where we discuss how the model reproduces the quantum 
predictions) . 

First, we will review the modifications to the equations that are caused by introducing 
the momentum cut-off. The most important difference is that the form of the guidance 
equations needs to be slightly modified in order to guarantee a pilot-wave dynamics 
which preserves the probability density p^. 

Then we will study what the typical beable configurations are corresponding to 
different macroscopic states. This is a prelude to the final subsection in which we study 
how our pilot-wave model reproduces the quantum predictions. 



4.1 Introducing the cut-off 

We introduce an ultra-violet momentum cut-off by making the replacement J d 3 p — > 
I\p\<a^ 3 P m ^ e Fourier expansion of the field operators. This is equivalent to putting 
the creation and annihilation operators to zero for momenta higher than the cut-off. As 
a result the field anti-commutation relations change to 

{Vv(x), $l, )0 ,(x')} = 5 AV cW5 (A) (x - x') , (42) 

where 

5( A )(x) = -i-/ dV p ' x - (43) 

( 27r ) J\p\<A 

Since the distribution j( A )(x) does not vanish for 0, the field operators ip and ip^ do 
not anti-commute anymore at different points in space. This causes some modifications 
in the equations of the pilot-wave theory, as we will soon see. 

First note that the distribution 5^ A ^(x) still acts as a 5-distribution on the space of 
functions for which the Fourier series only contain momenta |p| < A, i.e. if 

/(x) = - 7 =L=/ d 3 pf(p)e^, (44) 
V (2vr) 3 J\ p \<a 

then 

J rf 3 x/(x)5 (A) (x-y) = /(y). (45) 
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On the other hand, if we consider two such functions / and g which have a momentum 
cut-off, then their product fg generally has Fourier modes with momenta higher than 
the cut-off. This means that in general 

J d 3 */( x M x )<5 (A) ( x -y)± f(y)g(y) . (46) 

This property will become important when we search for possible guidance equations. 

By replacing the ordinary Dirac ^-distribution by the modified distribution S^ A \ most 
of the equations of Section 3 are immediately generalized to the case in which we have 
a momentum cut-off. One can still introduce the states |xi, Ai, a\, . . . ; x n , A re , a n ) by 
application of the operators i0' onto the state (O^), similarly as in (23). However, now 
the states are normalized with respect to the distribution 5^ A \ i.e. 

(xi, Ai, Ox! • • • i x ni A n , a n |x 1; A 1; a^; . . . ; x n , X n , a n ) 

= h. E(- 1 )^ (a) ( x -(d - oKvk^a ■ ■ ■ 5(A W) - 0Wv.)< • ( 4? ) 

Therefore different states are not orthogonal anymore (so that they form an overcomplete 
basis). As a result, the measure 

P(d 3 x 1 . . . d 3 x n ) = |xi,Ax,ai;...;x n ,A n ,a n ) ® |^) 

Ai,oi,...,An,a re ,| 

x (^| (g) (xi, A 1; ax; ... ; x n , A n , a n \ d 3 Xx . . . d 3 x n (48) 

is no longer projection valued but positive operator valued. 

Another consequence of the modified commutation relations of the field operators is 
that the states |x 1; Ax, ax; ... ; x n , A n , a n ) are no longer eigenstates of the fermion number 
density operator Fd(x) = J2x a a ( x )V ; A,a(x), although they are still eigenstates of the 

fermion number operator Fd, with eigenvalue n. 

Since the fermion number operator Fd commutes with the Hamiltonian, the fermion 
number is conserved. Assuming fermion number superselection as before, we can assume 
a Hilbert space which contains only states with a fixed fermion number. In the Hilbert 
space with a fixed fermion number n, a quantum state \tp(t)) can be expanded as 

1^(0) = ^2 J ■ ■ ■ d3xn 

Ai,ai,...,Anjttm£ 

x V'A 1) oi,...,An,on,f( x i) • • • > x n, t) |x 1; Ax, ax; ... ; x n , A n , a n ) ® |^) , (49) 

where the expansion coefficients 4)\ 1 ,a 1 ,...,\ n ,a n ,t('X-i, . . . , x n ,t) are given by 

V , Ai,oi,...,An,an,f( x i» • • • , x n , t) = ((xi, Ax, a x ; ... ; x n , A n , a n \ ® (f |) \ip(t)) , (50) 

just as in (28). From their definition it follows that the expansion coefficients are anti- 
symmetric under permutations of the indices 1, . . . ,n and contain only Fourier modes 
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with momenta smaller than the cut-off. The latter property can easily be seen to 
follow from the fact that the states |x 1; A 1; ai, ... ; x n , X n , a n ) only have Fourier modes 
with momenta smaller than the cut-off. This property will be important in evaluating 
integrals containing (^^-distributions. 

To construct a pilot-wave model we take the measure P(d 3 xi, . . . , d 3 x n ) as the po- 
sition POVM. For a state with a definite fermion number, the position density 
corresponding to the POVM P(d 3 Xi . . . d 3 x n ) is still of the form 

/(Xl,...,X n ,t) = ^ I^A 1 ,a 1 ,...,A„,a n ,c(xi,...,X n ,t)| 2 , (51) 

Ai,ai,...,A re ,an>£ 

where it is understood that the ?/>Ai,ai,...,A n ,a n ,£( x i, • • • > x n, t) now only have Fourier modes 
with momenta smaller than the cut-off, but the continuity equation for the density gets 

an extra term, which is due to the fact that the commutator P(d?x\ . . . d 3 x n ), Hj is no 

longer zero. The fact that the commutator is no longer zero follows from the property 
(46) above. 

In order to see what exactly happens, let us abandon the generality at this point 
and let us for simplicity just consider a single fermion. Similar expressions hold in the 
case of many fermions. For a one particle state the position density reads 

/(x,t) = ^ w (x,t)| 2 . (52) 

A,a,£ 

The position density satisfies the equation 



<V(x, *) 



dt 

cf. Section 2, from which it follows that 

dp* 
~dt 

where 



d 3 x+ (ip(t)\i P(d 6 x),H \i/j(t)} = 0, (53) 



V ■ (vV) + / = , (54) 



^ = E V>w a ~^WP* (55) 
is the usual velocity field and where 



g*(x,t) = l f d M^\K,a U A 2 ,a 2 (y)\t2)5 W (l 



Ai,Aa,ai,a2,Si,?a 

X W 1 ,ai,fi( x >*)V'A aia2l &(y,*)-V'A 1 ,ai4 1 (y»*)V'A 3l a ai f a (x,t)) , (56) 

is a contribution which comes from the term (ip(t) \ i[P(d 3 x), Hj] \ip(t)) in equation (53). 
This contribution does not vanish, since as mentioned above, the distribution <j( A )(x— y) 
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will not act as a 5-distribution, because we integrate over a product of two functions 
whose Fourier expansion has a cut-off. 7 

The term implies that the field can not serve as the velocity field for the beables 
anymore, since it would not preserve the density (r . One possible way to introduce a 
dynamics which preserves the density p^ is by supplementing the motion determined by 
the velocity field by a jump process, as explained in [23-27]. In this way the particle 
would move along a deterministic trajectory, determined by the velocity field v^ 1 , until a 
jump occurs. It is expected that the distance over which the particle jumps is typically 
of the order 1/A. 

Instead of introducing the jumps, which are stochastic, there probably also exist 
possibilities of keeping the dynamics deterministic. Suppose, for example, that g^ is 
such that we can write 

^(x,t) = V-(VV-V(x,t)) , (57) 

where 



J 47r|x-y| 

Actually checking that we can write g^ as in (57) would lead us too far, since g^ depends 
on the bosonic operators and hence on the details of the quantization procedure. If we 
can write g^ as in (57) then the equation (54) reduces to the continuity equation 

| + V.((v* + v*)/)=0, (59) 

where 

= V (V~V ) lp* . (60) 
Motivated by this continuity equation, we could introduce the guidance equation 



^ = v^ + v^. (61) 
at 

Since the contribution g^ disappears in the limit where the cut-off A is taken to 
infinity, we suspect that the importance of the contribution g^ for the dynamics of the 
beables decreases with an increasing cut-off A. Therefore we do not study the modified 
guidance equations in more detail. 

With the introduction of a cut-off the equations of motion for both the quantum 
state and the beables become cut-off dependent. The next step would be to remove 
the cut-off by applying a suitable renormalization scheme and by taking the cut-off to 
infinity. This would require a rather careful analysis, which is beyond the scope of 
this paper. One of the problems associated with taking the cut-off to infinity is that 
the number of position beables would go to infinity again (assuming that we have a 
bounded space), so that the problem of defining a measure on an infinite dimensional 
space reappears. On the other hand, there could be a natural ultra-violet cut-off, so 
that taking the cut-off to infinity would not be required. 



7 Note in passing that J d 3 xg^(x, t) = 0, so that the equation (54) implies that, after integration over 
the complete space, d(J d 3 xp^(x.,t))/dt = 0, which just expresses the conservation of fermion number. 
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4.2 Typical beable configurations for macroscopic states 

In order for our pilot- wave model to reproduce the predictions of standard quantum the- 
ory it is crucial that different macroscopic states correspond to different typical beable 
configurations. Having different typical beable configurations for different macroscopic 
situations namely implies that measurement results get recorded in the beable config- 
uration. In the following section, we discuss the equivalence of our pilot-wave model 
with standard quantum theory in some detail. In this section, we first consider the 
question what the typical beable configuration is for a macroscopic state and whether 
such typical configurations exist in the first place. 

In order to be precise, we say that a property of a beable configuration is typical 
if the property holds for almost all beable configurations. With 'almost all' we mean 
that the set of beable configurations which do not exhibit the property has very small 
p^(xi, . . . ,x n )<i 3 xi . . . G? 3 x n -measure, where is the position density. We say that a 
typical property holds for typical beable configurations. 8 

For example, in non-relativistic pilot- wave theory a typical beable configuration for a 
macroscopic system is such that the actual position distribution in physical space gives, 
on some coarse-grained level, the classical image of the macroscopic object. For example, 
a typical beable configuration for the state of a table (given by the wavefunction of the 
table) gives the classical image of a table. 

In the pilot- wave model presented in this paper we also have particles for macroscopic 
objects, like tables, but here these particles move against a background of 'vacuum 
particles' (of course, the distinction between particles belonging to the vacuum and the 
particles belonging to the system is only heuristic since the particles do not carry any 
label which would indicate such a distinction). As discussed in the previous section, 
this background of vacuum particles proves problematic if we do not introduce an ultra- 
violet momentum cut-off, since without cut-off the expected number of particles in any 
volume is infinite. It is only when we introduce a momentum cut-off that the expected 
number of particles becomes finite. 

Let us now study what a typical beable configuration will look like for a state de- 
scribing a macroscopic system. For the vacuum state, we expect that the typical beable 
configuration is such that it gives an approximate uniform distribution of positions 
in physical space. By 'approximate' we mean that the fluctuations in the distribution 
become negligible after suitable coarse-graining. For the state corresponding to a macro- 
scopic system, we expect a distribution which is made up of a uniform distribution of 
particles corresponding to the vacuum and a distribution of particles corresponding to 
the system, where the latter is localized around the location of the macroscopic system. 

Due to technical complications we do not rigorously prove that our expectation 
is correct. Instead we provide arguments that support our expectation. Let us first 
consider the vacuum state. The physical vacuum is not given by the state |0), which, 
in the language of standard quantum field theory, contains no particle or anti-particles, 
and which would be the vacuum if there were no interactions; due to the interactions 

8 The language of typicality is borrowed from Diirr et al. [38] who use it in the context of justifying 
quantum equilibrium. 
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the actual vacuum is much more complex. In the case there were only electromagnetic 
interactions, the physical vacuum could be thought of as some superposition of the 
state |0) with states containing one particle-anti-particle pair, states containing two 
pairs, and so on. However, it is unclear to us what exactly the physical vacuum looks 
like if also the weak and strong interactions are brought into account. Therefore instead 
of considering the physical vacuum, we consider the state |0) for simplicity. For the 
state |0) the expected number of fermions in a volume B, which we denote by uq(B), is 
given by (see Appendix B for the details of the calculation) 

n (B) = (0| F d (B) |0> = £ / d 3 p [ Sx = ^A 3 V(B) , (62) 

where V(B) is the volume of the region B. This means that the expected fermion 
number density is uniform. In order to find out whether this behaviour is typical for 
beable configurations corresponding to the state |0), we calculate the standard deviation 

A F d (B) = J <0| F d {Bf |0> - (0| F d (B) |0> 2 = y/ <0| F d {Bf |0> - n (B) 2 . (63) 

The standard deviation gives an estimate of the possible fluctuation of the fermion 
number in the region B. A small standard deviation would indicate that a uniform 
distribution is typical. 

In order to calculate the standard deviation, we assume a spherical region with radius 
b, which we denote B s . We consider radii b which are of the order of the atomic distance 
10~ 10 m = 1A or larger. This means that b is at least two orders of magnitude larger 
than the Compton wavelength of the electron, which is of the order 10 _12 m. Hence b is 
also much bigger than the Compton wavelengths of the quarks and the other electrically 
charged leptons. At this point, we also assume the cut-off to be of the order of the 
Planck mass, i.e. 1.3 ■ 10 19 GeV ~ 10 35 /m. 

The details of this calculation are given in Appendix B. The result is that the 
fluctuation of the number of fermions in B s is 

A F d (B s ) « 0.39v / AV(5 s ) 1 / 6 . (64) 

If the fluctuations in different regions were uncorrelated then according to classical 
statistical physics one would expect the fluctuation to be of the order of n (B s Y^ 2 « 
0.llA 3 / 2 V(B s y/ 2 , i.e. the square root of the number of particles in the region B s [46, 
pp. 8-9]. Since A F d (B s ) « 0.39^/AV(B s ) 1 ^ 6 the fluctuation per unit volume decreases 
faster with increase of the volume than expected in the absence of correlations. So the 
fluctuations must be correlated. 9 

Since we do not consider the physical vacuum, but only the free vacuum |0), the only 
cause for suppression of fluctuations comes from the mutual repulsion of the vacuum 

9 Similar conclusions were drawn by Hciscnbcrg [47]. Hciscnbcrg performed a calculation similar to 
ours. The main difference is that Hciscnbcrg also introduced a coarse-graining over time (which signif- 
icantly simplified the computations). He had the additional result that the fluctuation also decreases 
when there is a coarse-graining over larger time intervals. 
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particles due to the Pauli exclusion principle. It is unclear what exactly would happen 
in the physical vacuum. If there were only electromagnetic interactions and only one 
type of charged particle then the fluctuations would be even more suppressed because of 
vacuum polarization effects. However, it is unclear what exactly happens if we include 
all types of particles and interactions. 

Before giving a numerical estimate for the fluctuations, we first consider a state 
describing a macroscopic system. Let us assume a state of the form 



W) = -7=f Y] ( d3 Pl ■ ■ f d i p m ip Xl ,s 1 ,...,X m ,s m (Pl, ■ ■ ■ , Pm) 

X u s u ...,X m ,s m J \Pi\<A J\Pm\<A 

xA 1 ,sM---A m ,sJPm)\0) , (65) 

where the expansion coefficients </?Ai,si,...,A m ,s m (pi, • • • , Pm) are given by 

<f\ usl ,...,x m ,s m (-pi,---,Pm) = -7==(0|c Aro , Sm (p m ) . . .c Xl , sl (pi)\cp) (66) 



when all momenta are smaller than A and by zero otherwise. Note that the state contains 
m particles on top of the vacuum |0) (not the vacuum |0u)) and no anti-particles. We 
will further assume that the state contains only 'ordinary matter', i.e. up and down 
quarks and electrons (in fact we should also include the background of neutrinos, but 
as we will argue below, the particle number density of neutrinos is considerably smaller 
than the particle number density of the ordinary matter, so that this background can 
be ignored for our purposes). 
Using the Fourier expansion 



V 2 Ai,ai,...,A m ,a m ( x l> • • • 5 x m) 

d 3 pi f d 3 p m I m Xl 



Ipi|<a yWf"'J\ Pm \<A x/l^fy E Xl (p 



-u Xl , Sl>ai ( Pl )e^ 



J 1 )™ , ux m , Sm ,aJp m )e iPm - Xm ip Xl , Sl ,..., Xm , Sm (pi, ■ ■ ■ , p m ) , (67) 

the state can be written as 

\ { P) = -r=x I d 3 x 1 ...d 3 x m (p Xltair _., Xmjam (x 1 ,...,x m ) 

Vffl! x — ; J 

Ai,ai,...,A m ,a m 

x^ 1 ,a 1 (xi)...^ m , am (x m )|0) , (68) 

where 

^Ai,ai,...,A m ,o m ( x l5 • • • ) x m) ,a m ( x m )...^A l! ax(xi)b>. (69) 
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Assuming the state \cp) to be normalized to 1, we have 



,A miSm (Pl,---,Pm)| 2 - (70) 

As with the vacuum state, our representation of a macroscopic system by the state \tp) 
is not entirely accurate. For the state \<p), the only interaction between the macroscopic 
system and the vacuum is caused by the Pauli exclusion principle. In general, we also 
have to include the standard interactions between the system and the vacuum, like the 
electromagnetic interaction, which will cause a reorganization of the vacuum particles. 
Again we only know what would happen if we only had electromagnetic interaction. In 
that case, we would have, in the standard particle-anti-particle picture, the appearance 
of virtual particle-anti-particle pairs around electrically charged particles. 

We want to calculate the expectation value and standard deviation for the fermion 
number for a region B of the order of the atomic distance 10 _10 m = lA or larger, just as 
we did for the vacuum |0). In order to simplify the analysis, we assume that the state \tp) 
represents particles that are approximately localized within this region B. This means 
that most of the support of the wavefunction f\ ltai ,...,x m ,a m i^-i, ■ ■ ■ , x m ) is within the 
region B m C IR 3m . Note that the wavefunction cannot have all of its support within the 
region B m , because it is a superposition of positive energy states (i.e. a superposition of 
the functions u(p)). At best, most of the support can be within a region of the order of 
the Compton wavelength of the lightest particle under consideration, which in this case 
is the electron Compton wavelength 10 _12 m (we are assuming that the state represents 
ordinary matter). The effect of the cut-off on the localizability of the wavefunction is 
hereby negligible, since it only becomes important on scales of the order of 1/A, which 
is the Planck length when A is the Planck mass. Since we are interested in regions B 
with a radius of the order of the atomic distance 10 _10 m = lA or larger, the fact that 
most of the support can at best be within a region of the order of the electron Compton 
wavelength, will only be of minor importance. 

The details of the calculation are given in Appendix C. We have that the expected 
number of particles in the region B is given by 

(<p\F d (B)\ip)&m + n (B). (71) 

For a region B C 1R 3 \ B we have that 

(cp\ F d (B) \cp) « n (B) . (72) 

So in the region B the expected number of particles is approximately given by m, the 
number of particles corresponding to the macroscopic object, plus uq(B), the number 
of vacuum particles in the region B. In a region B outside B the expected number of 
particles is given by the number of vacuum particles in that region. 



1= Yl d 3 Xl---d 3 X m \(p Xl ,a 1 ,...,X m ,a m (^U 



Ai,ai,...,A m 

£ 

Ai ,si ,...,A m ,s m 



d pi ■ ■ ■ / d p m \{p\ 1)S1 , 

pi|<A ^|Pm|<A 
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The standard deviation is given by 



A v F d (B) = ^{ v \F d {BY\ v )-{ v \F d {B)\ v f » A Q F d (B) . (73) 

For a region B outside B, we have the same expression for the standard deviation. 

Let us now consider the question whether states corresponding to macroscopically 
distinct systems will correspond to different typical beable configurations. Two states, 
say and |</> 2 ), will correspond to distinct beable configurations if there exists a region 
B for which 

| (</?i| F d (B) |^) - (ip 2 \ F d (B) \<p 2 ) | > m a x(A Vl F d (B), A V2 F d (B)) « A F d (B) . (74) 

In order to address the question let us consider a specific example. Consider two states, 
one which represents solid matter, say graphite, in a spherical region B s with radius b, 
and one which represents, say, air in that region. What should now be the volume V(B S ) 
of the region B s in order for the states to correspond to different typical fermion number 
configurations in B s l Answering this question will give us an idea whether macroscop- 
ically distinct systems will correspond to different typical beable configurations. If we 
find that we only need a small volume V(B S ), e.g. microscopically small, then this would 
suggest that the answer to the latter question is affirmative. Now, since the state of 
the solid matter has a considerably higher fermion number density in B s than the state 
representing air, the volume V(B S ) should be such that 

m > A F d (B s ) , (75) 

where m is the number of fermions corresponding to the solid matter contained in 
B s . The fermion number density of graphite approximately given by 4.2 • 10 30 /m 3 (1 
carbon atom consists of 42 fermions, namely 18 up quarks, 18 down quarks and 6 
electrons; 0.012kg carbon contains 6 ■ 10 23 atoms and the mass density of graphite is 
approximately given by 2-10 3 kg/m 3 ), so that m ~ V(B S )4.2-10 30 /m 3 . We hereby ignored 
the background of neutrinos, since the density of neutrino is considerably smaller than 
the fermion density of graphite. 10 Taking the cut-off A of the order of the Planck mass 
10 35 /m, we have that (75) is valid when the radius b of B s satisfies b 3> 2.6 ■ 10~ 6 m. This 
is a rather large lower bound (for comparison the diameter of a human hair typically 
ranges from 10 _4 m to 10~ 5 m). However as mentioned before, this result would be 
modified when all the interactions would be taken into account. If there were only 

10 The neutrino density on earth is believed to be dominated by relic neutrinos, i.e. the neutrinos 
which decoupled within the first 10 seconds after the Big Bang. Cosmological considerations lead to 
an estimate of a neutrino density of approximately 3.4 • f0 8 /m 3 [48, p. 222]. For comparison, the sun, 
which is the main steady source of neutrinos for the earth, gives a density of approximately 3.3 • 10 6 /m 3 
[49, p. 14]. The rest of the stars of our galaxy together only yield a density which is at most 10 -8 
times the density of solar neutrinos [49, pp. 165-166]. Other sources of neutrinos are supernovae. The 
well-known supernova SN1987A in the Large Magellanic Cloud, a nearby dwarf galaxy, is estimated to 
have produced a density of at most 10/m 3 , and this over a period of less than 15 seconds [49, pp. 443]. 
A typical supernova in our galaxy is estimated to produce a density which is 25 times larger [49, pp. 
437]. 
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electromagnetic interactions and one type of particle then we know that the result 
would improve. 

Note also that the estimate is cut-off dependent. A lower cut-off would lead to a 
decrease of the lower bound. A higher cut-off on the other hand would lead to an 
increase of the lower bound. This seems to imply that taking the cut-off to infinity in a 
renormalization process would obscure the identification of macroscopic objects in the 
beable configuration. On the other hand, in the case there is a natural cut-off, which is 
not too high, there is no problem with such an identification. 

4.3 Reproducing the quantum predictions 

Given that macroscopically distinct states imply distinct typical beable configurations, 
it is straightforward to show that our pilot-wave model reproduces the predictions of 
standard quantum field theory. The arguments are similar to the ones used to show 
that non-relativistic pilot-wave theory reproduces the predictions of standard quantum 
mechanics. For definiteness we repeat these arguments here. We closely follow the 
original exposition of Bohm, who was the first to present a convincing analysis for non- 
relativistic pilot-wave theory [2] . 

Since quantum theory makes probabilistic statements about outcomes of measure- 
ments, we just consider the measurement situation in detail. According to standard 
quantum theory, or at least according to some versions of it, a measurement situation 
can be described as follows. The wavefunction ip (which is the expansion coefficient as 
defined in equation (49)), which describes the system under observation together with 
the measurement device, macroscopic pointer and environment, evolves into a superpo- 
sition X^i c iV'i) where the ipi are normalized to f. The different ipi correspond to the 
different possible outcomes of the measurements, i.e. the different ipi describe the dif- 
ferent possible macroscopic pointer configurations which indicate the different possible 
measurement outcomes, like for example a needle pointing in different directions. The 
wavefunction then collapses to one of the in the superposition, say ipk, and this with 
probability |cfc| 2 . The state ipk represents one of the possible macroscopic pointer config- 
urations and hence also the corresponding outcome of the measurement. So the collapse 
ensures that measurements have definite outcomes, whereby the outcome is determined 
by the wavefunction through the eigenstate-eigenvalue link. Of course, as is well known, 
it is rather ambiguous when the collapse exactly occurs and this ambiguity forms the 
core of the measurement problem. 

In our pilot-wave model the wavefunction evolves according to the Schrodinger equa- 
tion at all times. The wavefunction never undergoes collapse. Outcomes of measure- 
ments are not recorded in the wavefunction, like in standard quantum theory, but in 
the beable configuration. This is seen as follows. In a measurement situation the wave- 
function evolves into a superposition ^jQ^, just as in standard quantum theory, but 
there is no subsequent collapse. Since the ipi correspond to different macroscopic situa- 
tions, like for example a needle pointing in different directions, they will have different 
typical beable configurations. This implies that the vpi have minimal overlap, i.e. that 
they have approximately disjoint supports in configuration space. This means that the 
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actual beable configuration (x 1; . . . , x n ) will be in the support of only one of the wave- 
functions ipi, say tpk, or in other words the actual beable configuration will only be a 
typical configuration for one of the wavefunctions In this way, the beable configu- 
ration displays the measurement outcome. Since in equilibrium, the probability for the 
beable configuration (xi, . . . , x n ) to be in the support of tpk is given by |cfc| 2 , we recover 
the quantum probabilities. 

Further, if the different wavefunctions ipi stay approximately non-overlapping at later 
times, then one can easily check that only the wavefunction ip^ determines the velocity 
field for the beable configuration. Since the other packets ipi, with i ^ k, play no 
role in determining the velocity field anymore, one can simply remove them from the 
description of the beable configuration. This is called an effective collapse and clearly 
explains the success of the ordinary collapse in standard quantum theory. 

5 Note on the predicted violation of fermion num- 
ber conservation in the Standard Model 

The assumption that fermion number F is conserved is a key ingredient that allows us 
to construct our pilot-wave model. This assumption is justified in the energy regime 
where the electroweak gauge symmetry is broken. In the same regime, one actually has 
that both quark number Q and lepton number L are conserved. The quark number is 
defined as the number of quarks minus the number of anti-quarks Q = n q — rig. The 
lepton number L = rii — n\ is the number of leptons minus the number of anti-leptons. 
The fermion number F is the sum of the two, i.e. F = Q + L. Conservation of both 
quark number Q and lepton number L actually allows us to also introduce additional 
beables corresponding to the particle type, i.e. quark or lepton. In such a pilot-wave 
model the actual particles would not only have a position but also a label indicating 
whether it belongs to the family of quarks or the family of leptons. 

However, the Standard Model predicts a violation of conservation of all three num- 
bers F, Q and L for sufficiently high energies. The violations of conservation of these 
quantum numbers are so-called anomalies. 11 This means that classically, i.e. before 
quantization, these quantum numbers are conserved. It is only after quantization that 
they are not conserved anymore. Before quantization the theory is invariant under the 
global U(l) symmetries ip\ —>■ e ta ip\, where the index A either ranges over all possible 
values, or only the values corresponding to quarks, or only the values correspond to 
leptons. By Noether's theorem this respectively gives us conservation of F, Q and L. 
When quantizing one needs to give up either these global symmetries or the electroweak 
gauge SU(2) x U(l)y symmetry. More specifically, it is the regularization that goes 
with the quantization that requires that either of these symmetries is given up. Since 
the gauge invariance is important to show renormalizability, the global U(l) symmetries 
are given up, so that the quantum numbers F, Q and L are no longer conserved. 12 The 

11 A clear introduction to anomalies is given in [50], [51, pp. 1-102] and [52, pp. 217-235]. 

12 For example, if we were to apply a similar cut-off regularization as in Section 4 to the 'full-blown' 
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possibility of non-conservation of these quantum numbers is important in cosmology 
since it might yield an explanation for baryogenesis, see e.g. [53-55]. 

Interestingly, often a heuristic derivation of anomalies is given in the Dirac sea picture 
(see e.g. Jackiw [50, pp. 333-337], [51, pp. 22-26], [56], and Greiner et al. [52, pp. 219- 
221]). Maybe this heuristic derivation sheds a light on how to construct a pilot- wave 
model in the case when fermion number conservation is violated. 

An alternative would be to consider the quantum number B — L where B is the 
baryon number B = Q/3. This quantum number is conserved in the standard model. 
It is also conserved in unification theories like e.g. the 57/(5) unification theory [57]. 
Theories that conserve B — L also go by the name B — L theories [57]. Maybe one can 
introduce beables by exploiting the conservation of the quantum number B — L. 

6 Outlook 

We further developed the Dirac sea pilot-wave model originally presented by Colin and 
Bohm et al. Several points need further attention. First of all, the model relies on the 
conservation of the fermion number F. Since the Standard Model predicts a violation of 
fermion number conservation for sufficiently high energies, our pilot-wave model should 
be further developed in order to be able to deal with this possibility. It is unclear how 
this could be done in a way that preserves determinism. An alternative would be to try 
to introduce beables for another quantity which is conserved according to the Standard 
Model, like baryon number minus lepton number B — L. 

We have further shown how the pilot-wave model can be regularized by introducing 
an ultra-violet momentum cut-off and by assuming a finite space. In this context we 
also discussed the emergence of the familiar macroscopic world in the model. Although 
we gave arguments to support our expectation that the familiar macroscopic world does 
emerge, a more conclusive analysis would be desired. 

In our model no beables were introduced for the bosonic degrees of freedom of the 
quantum state. If the beables in our model provide an image of the familiar classical 
world, as we expect, then there is actually no need to introduce beables also for the 
bosonic quantum fields. Bell, for example, did not introduce beables for the bosonic 
degrees of freedom in his lattice model. On the other hand, one could introduce beables 
for the bosonic degrees of freedom by using the field beable approach initiated by Bohm 
for the quantized electromagnetic field. Another possibility would be to introduce parti- 
cle beables also for the bosons, although it is unclear for the moment how this could be 
done in a satisfactory way. In connection to this, it might also be interesting to look at 
the supersymmetric version of the Standard Model. Since the supersymmetry implies 
a symmetry between bosons and fermions, maybe an ontology which treats bosons and 
fermions in a similar way suggests itself in this context. 

electroweak theory (i.e. not merely the spontaneously broken sector of the electroweak theory), then 
the fermion number could still be conserved since this type of regularization breaks the electroweak 
gauge symmetry 
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A Commutator of the position PVM and the inter- 
action Hamiltonian 

We show that P{d 3 x\ . . . d 3 x n ) commutes with Hi, where P(d 3 x\ . . . d 3 x n ) and Hi are 
respectively given in (27) and (13) (no cut-off is introduced here). 

In order to do so, first note that if an operator A commutes with Hi and with 
h\. a ,\' ,a'( x ) for all values of A, A', a, a', x, then a a ( x )^' ? /'A i a( x ) commutes with Hj. 
This follows easily by using the commutation relations 

^A,a(x),5j] = ^/iA,a,A',a'( x )^A',a'( x ) , 
A', a' 

A', a' 

If we now take A = \0d) (g> |£) (£| <g> (0d\, then A commutes with Hi and with 
^A,a,A',a'( x ) fo r an values of A, A', a, a', x so that JTJ A a ip\ a (xi)Aipx,a(^i)d 3 xi = P(d 3 x\) 
commutes with Hj. Then by induction and the fact that 

P(d 3 X 1 . . . d 3 X n ) = ~J2 ^A.aM^^l ■ ■ ■ d 3 Xn-^x^^d 3 X n , (77) 

A, a 

we find that P(d 3 x± . . . d 3 x n ) commutes with Hi. 

As we explain in Section 4, this property no longer holds if a momentum cut-off is 
introduced. The reason is that (76) no longer holds in that case. 

B Expectation value and standard deviation of the 
fermion number for the vacuum |0) 

In this appendix, we calculate the expectation value 

n (B) = (0\F d (B)\0) = [ d 3 x(0|F d (x)|0) (78) 

JB 
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and the standard deviation 

A F d (B) = yj (0\F*(B)\0) - MB) 2 (79) 

for the fermion number in a region B for the vacuum |0). In the calculations we use 
an ultra-violet momentum cut-off |p| < A. The cut-off is assumed to be at the Planck 
scale. 

The expectation value for the fermion number in the region B is given by 



A.o 



(B) = J2 d 3 *(0|^ a (x)^, Q (x)|0> 
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P\3 J2 I d3x I d3 Pi I d3 P2- 
ln ) , TT 1 . JB JlniKA i| P2 l<A 



^ X,a^,sJ B J \^ A ^IP2l<^ \/^a(Pi)^a(P2) 

(0MA, Sl (pi)rfl, S2 (p 2 )|0)< Sl , a (pi)^ 2 , a (P2)e^- P2 )- x 



A,.s 

A3 £ / - 



3tt 2 

A 



P|<A 

3, 



X 

B 



= -A 3 V(B) , (80) 
where we have used the normalization of the spinors 

J2 V {, S1 ,a(P) V \s2,a(P) = ^s 1 s 2 Ex(p)/m X (81) 
a 

in the third equality. In the last equality we used the fact that Yl\ 1 = 24 (we have six 
types of quarks, which can each have three different colours, and six types of leptons). 
We also introduced the notation V(B) for the volume of B. Since the expectation value 
is proportional to the volume of the region B, the expected fermion number density is 
uniform. 

In order to get the standard deviation, we first calculate 

(0\F!(B)\0)= f d * x I A(0|^, a (x)^A,a(x)^, a ,(y)^,a'(y)|0). (82) 

A,A' ,a,a> J B Jb 

Since any field operator ip is a superposition of operators c and eft, there will be two 
non-zero contributions to this expression: one coming from the product deft deft and one 
coming from the product We call these contributions respectively Cj and Cu, 
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so that (0|F|(S)|0) = d + C J7 . We have that 

^ = TTT^TTT ^ I d 3 x I d 3 y [ d 3 p 1 d 3 p 2 d 3 p 3 d 3 p± 
\ 2n ) x fri, . Jb Jb J\ Pk \<A 



V / j Ba(pi)^a(p 2 )^v(P3)^v(P4) 
Using the canonical anti-commutation relations (17), we find 

(0|rfA, Sl (Pl)4, S2 (P2) rf A',s 3 (P3)4', S4 (P4)|0) = 5 ain 5 SiSA 5(pi - p 2 )5(P3 - P4) , 

so that after integrating over p 2 and p 4 , we get 

Cl = lh« E / d ^ / ^ / / d 3 y = n (B) 2 . 
y 27T ) \&^.J\pi\<a J\ P3 \<A Jb Jb 



A, A', si, S3 



Now we compute Cu- We find that 



C u = TTTTfi 5^ I d 3 x I d 3 y [ d 3 p l d 3 p 2 d 3 p 3 d 3 p A 

-7=7 =^S= (0Ma,s 1 (Pi)ca ] s 2 (P2)4,s 3 (P3)^,s 4 (P4)1Q) 
V/^A (pi ) £ A (p 2 ) £y (p 3 ) £? A / (p 4 ) 

X^,s 1 ,a(Pl)«A,s 2 ,a(p 2 )4,S3,a'(P3)^,s 4 x(p4)e^ Pl+P2 > X e-^ + ^)^ 

Using 

(0Ma,si(Pi)cA !S2 (P2)c A ', S3 (P3)4',S4(P4)|0) = ^AA'5s 1 s4^ 2 s 3 5(Pl ~ P4)<$(p 2 ~ Ps) 

and integrating over p 3 and p 4 , we get 

CU = a JL L *° Sb ** L, L, 

<s 1 ,a(Pl)^ !S2 ,a(p 2 )4 S2 ,a'(P2)^ Sl ,a'(Pl)e i(Pl+P2Hx - y) 



= tAg J2 I d3x I d3 y I ^ I rf3 ^ 2 

[2'k) [0 ^J b J b J\ Vi \<a J\ P2 \<a 

E X (pi)E X (p 2 ) + pi ■ p 2 - 7»A i(p 1+P2 )-(x-y) 

Ea( Pi )^a(p 2 ) 
In the last equality, we have used the relations 



/ x t / \ ^a(p) +a-p + /9m A \ 

E"WpR,, a -(P) = ( gA(p)+ 2 a mA p " /3mA ) , 
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So finally we obtain (0|F|(B)|0) = d + C n = n (B) 2 + C n , so that (A F d (B)) 2 = 

(0\F»{B)\0) ~ MB) 2 = Cn. 

We now evaluate this expression for a spherical region B, which we now denote by 
B s , with radius 6 and centered around z. We are interested in radii 6 which are of the 
order of the atomic distance 10~ 10 m = lA or larger. 

In order to simplify the equations we will make the replacement 

o / 2 f n |X-Z| 2 

fc/ d 3 xe~^" , (90) 

where the factor a/2/97t is a normalization constant. In this way we get, after Gaussian 
integration over x and y, and after change of integration variables (pi,p2) — * (p, q) 
where Pi = p and P2 = — p — q, that 1 



.13 



(AoWf d3 P [ 



£a(p)£a(p + q) - P • (p + q) - m\ 

^A(p)^A(p + q) • 1 J 

We wrote an approximation sign instead of an equality sign because we dropped the 
integration bound |q + p| < A for q. This can be done because of the following reason. 
The main contribution to the integral comes from momenta q for which |q| < 1/6, hence 
the bound |q+ p| < A only plays a role when |p| < A — 1/6. Since A and since 

the integration domain of p is given by |p| < A, the removal of the bound |q + p| < A 
for q in the integral forms a good approximation. (Actually, more precisely, the right 
hand side forms of (91) forms an upper bound for (AqF^Bs)) 2 , since the integrand is 
strictly positive, as can easily be checked.) We can further integrate this by passing to 
spherical coordinates, (p) — > (p,(j> p ,9 p ) and (q) — > (q,(f) q ,9 q ), where we choose 9 q to be 
the angle between p and q. We obtain 

26 6 r+°° 



x 

oz,6 /"A r+oo 



(P,q), (92) 



where 



hip, q) = PQ (pq - \Ex(p + q)E x (p) + (g 2 - M)) (93) 

and E\{p) = y/p 2 + m\. Note that (A F d (B s )) 2 is independent of the centre z of 
the spherical region B s , which was to be expected since the vacuum |0) is translation 
invariant. 



13 The symbol is used for 'approximately equal as'. The symbol '<~' will be used for 'of the same 
order as'. 
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We see that (A F d (B s )) 2 is the sum of contributions of the different types of particles, 

(A F d (B s )) 2 = J2^oF d (B s ))l, (94) 

A 

where 

(A F d (B s )) 2 x = — J dp J dqe- b2 « 2 I x (p, q) . (95) 

In order to evaluate the integrals in (A F d (B s )) 2 x , we consider the following separate 
cases: 

• Case 1: b^> l/m\ 

• Case 2: b < l/m\ 

• Case 3: Otherwise, i.e. (b > l/m A ) A (-1(6 > l/m A )) 

Remember that we are interested in radii b which are of the order of the atomic distance 
10 _10 m = lA or larger. Therefore, the charged fermions, i.e. all types of fermions except 
neutrinos, always fall under Case 1. This is because the charged fermion with the largest 
Compton wavelength l/m\ is the electron, for which l/m\ ~ 10 _12 m. Cases 2 and 3 
apply to neutrinos. 

Case 1: b ^> l/m\ 

The main contribution to the integral comes from momenta \q\ for which \q\ < 1/6. Since 
we also assume that the radius b of the region B s is much bigger than the Compton 
wavelength, we have that the main contribution to the integral in (AoF d (B s )) x comes 
from momenta \q\ for which \q\ <g m\. By performing a Taylor expansion of E\(p + q) 
around q = 0, we find up to third order in \q\/m\ that 

hip, q) « ~Ex(p) 2 pq + ^(2/? 4 + 3p 2 q A m 2 x ) . (96) 

The first term will not contribute to the integral in (A F d (B s )) 2 x A since it is odd in q. 
Doing the integration over the second term yields 

(AoF^B^l^^-^Ab. (97) 
In this expression, we have dropped the terms of the order bm\ and lower since rri\ A. 

Case 2: b < l/m\ 

In order to evaluate (A F d (B s )) 2 x this time we write 

(A F d (B s )) 2 x = (A F d (B s ))l a + (A F d (B s )) 2 XA , (98) 
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where 

26 6 



26 f a r+°° 
(A F d (B s ))l a = — dp dqe- b ^I x {p,q) , 

J7r Jo J-oo 
niS /"A r+oo 

(A F d (B s ))l A = / dp dqe^h^q). (99) 



We take a such that A ^> a ^> 1/6. Together with 6 < 1/tha we have that a ^> 1/b > 
m x . 

Let us first evaluate (A F d (B s )) 2 x A . Since the main contribution to the integrals 
comes from momenta |g| for which \q\ < 1/6, the particular choice for a implies that 
the main contribution to (A F d (B s )) 2 x A comes from momenta p for which p ^> m x and 
[p + <z| 3> m^- For p ^> rrtA, |p + q\ 3> and (p, g) G [0, a] x [—1/6, 1/6], we have that 

h(p,q)~lq*-lp 3 q. (100) 

The second term will not contribute to the integral in (A F d (B s )) 2 x A since it is odd in 
q. Doing the integration over the first term yields 

(A F d (E s ))^ T ^A6. (101) 

This is the same expression as in (97). 

We now show that (A Fd(B s )) 2 x a <C (AoF d (B s )) 2 x A . In order to do so, first note that 
for (p,q) G [0, a] x [—1/6, 1/6] we have that 

\h(p, q)\ < PV + \p\q\E X {p + q)E x {p) + ^fc^ {\<1?P + H V) 

2 1 
< pV + ^P\q\E x (p + 9)^a(p) + g^Ab + q) {\q\ 3 + |g| 2 p) 



1 / / \\ 2 2 { 1 a 



3T 



(102) 



where in the last equation we have used the fact that a ^> m x and a> 1/6. As a result 
we have 

(Ao^(5 s ))L<i« 4 ^- (103) 



Now we check that (A F d (B s )) 2 x a is negligible compared to (Ao F d (B S )) 2 X A . To this end 
assume a ~ 10 3 /6 (a was chosen such that A ^> a ^> 1/6 > t^a). We then have that 
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(A F d (B s )) 2 Xa < 2 ■ 10 12 /9tt 3 / 2 ~ 10 10 . If we assume the cut-off to be of the order 

of the Planck mass, i.e. A ~ 10 35 /m, then (AoF d (B s )) 2 Xa is clearly much smaller than 

(A F d (B s ))l A « 6A/18tt 3 / 2 ~ 10~ 2 A6 ~ 10 33 6/m for b > 10~ 10 m = lA. Hence we 
obtain 

(A F d (B s )f x « (A F d (B s ))l A « y^A6 . (104) 

Case 3: Otherwise, i.e. (6 > l/m x ) A (-1(6 ^> l/m\)) 

This case corresponds to b larger than the Compton wavelength of the fermion but not 
too much larger either. Therefore we take b = 10 13 /rri\, where we take j3 £ [0, 3]. 

In this case, we use the same technique as in the previous case. We split up the 
integral like in equation (98). We choose A ^> a ^> m\. Together with the conditions on 
b, we have then that a ^> m\ > 1/6 > I0~ 3 m\. Just as in the previous case, it follows 
that the main contribution to the integral in (A F d (B s )) 2 x A comes from momenta p 
for which p ^> m\ and \p + q\ ^> m\. Therefore we have the same expression for 
(A F d (B s )Y XA and the same bound for (A F d (B s )) Xa . 

We also have that (A F d (B s )) 2 x is dominated by (A F d (B s )) 2 x A . To see this, consider 
the ratio e = (A F d (B s )) X a / (A F d (B s )) 2 X A w 4a 4 6 3 /A. If e < 1 then (A F d (B s )) 2 x 
will be dominated by (A Q F d (B s )) 2 x A . Now, since A ^> a ^> m\ > 1/6, we can take 
a ~ 10 3 mA. In this way, we get e < 4 • 10~ 23+3 ^mAm. Since the Case 3 only applies 
to neutrinos, we can find an upper bound for e by substituting m\ by the upper bound 
for the neutrino mass. Cosmological considerations imply that the sum of the masses 
of the different neutrino flavours is smaller than 2eV ~ 4 • 10 6 /m [32, p. 474]. Therefore 
an upper bound for e is given by 1.6 • 10~ 16+3/3 . Clearly this upper bound is still much 
smaller than 1 for (3 £ [0, 3]. Therefore we have, just as in the previous cases, that 

(A,F d {B s )) 2 XA ^^j- 2 Kb. (105) 

Summary 

Combining the above results we find that 

1 4 /4A 2 / 3 1 

(A F d (B s )) 2 = £ (A F d (B s )) 2 x « £ ^Ab = ^Ab = (-) ^AV{B S )^ , 

A A 

(106) 

where V(B S ) = 4-7r6 3 /3 is the volume of the spherical region B s . The standard deviation 
now reads 

A F d (B s ) « f^j /3 -l^A^ViB^ « 0.39A 1 / 2 V(5 S ) 1 / 6 . (107) 
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C Expectation value and standard deviation of the 
fermion number for a macroscopic state \ip) 

In this appendix, we calculate the expectation value and the standard deviation for the 
fermion number in a region B which is of the order of the atomic distance 10~ 10 m = lA 
or larger, for the state \ip) given in (65), which represents m particles (and no anti- 
particles) on top of the vacuum |0) that are approximately localized within B. Before 
we calculate the expectation value and the standard deviation, we first introduce some 
notation and useful identities, as well as some approximations. 



C.l Some notation and useful identities 

We first introduce the operator 



=T l rf3a; l--- rf3a; ™V 3 A 1 ,a 1 ,...,A m ,a m (xi,...,X m ) 

ml \ — : J 

i,ai,...,A m ,a m 

x ^! 1 «( x l)•••^'L,«„( , 

- y 



Ai,si,...,Ak 



d 3 p 1 ... / d 3 p m ip Xl ,8 1 ,...,\ m ,8 m (Pl, ,Pr 

Pi|<A J\p m \<A 



xcl llSl (Pi)---cL 



[Pr, 



for which Cj|0) = \tp). 
We have that 

where 



F d (B),dt =mdUB) 



'1081 



(109) 



Cl {B) = — = Y] / d 3 x 1 / d 3 x 2 . . . d 3 x m <p Xl , ai ,...,x m ,a m (xi, . . . , x m ) 
\jm\ . , Jb J 

Ai,ai,...,A m ,a m 

X V5 1 ,a 1 ( X l)---fL,a ro ( X m) 

This follows from 

[F d (B)^{ i ^ 1 )...^ a J^ m )_ 

m „ 

= / d 3 x5^(x - *)$ liG1 (xi) . . . ^U(x,) -> ^(x)] . ..^W 



(110) 



i=i 



'111] 



where 



We further have that 



means that ip\. .(xj) is replaced by .(x). 



F d (B),CUB) =mCJB), 



(112) 
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where 



CJB) = — i= [ d 3 X 1 [d 3 X 2 ...d 3 X m (f Xl) a 1) ...,X m ,a m (^l^--^m) 



X 



Ai,ai,...,A m ,,a 



??). 



-^^(xi) + ^ / fi ^I ijai (x)5^(x - Xl )) « 2 , a2 (x 2 ) . . .^L )am (x m ) . 

(113) 



Other useful relations are 

a v iv>> = a^at |o) = io) , 

C v Cl(B)\0) = (015^(5)10)10) = (<p\dl(B)\0)\0) . (114) 

The first relation can be derived as follows. The state C v \(p) = C^C^ |0) contains the 
product c. . . cc" . . .c 1 of m annihilation operators and m creation operators acting on 
the state |0). Moving the annihilation operators to the right, we find that C v \tp) ~ |0) 
and hence, by using the resolution of the identity, we have C v \<p) = |0)(0|C V ,|</?) = 
\0){ip\ip) = |0). The second relation can be derived completely analogously 

C.2 Some approximations 

We now present some approximations which we will use in the calculation of the ex- 
pectation value and the standard deviation. The approximations have to do with the 
fact that the wavefunction </?Ai,ai,...,A m ,a m ( x i ; • • • ; x m) only has most, and not all, of its 
support in the region B (as explained in Section 4.2). 
First, we have 

(V\CI{ B )\V>) = ~7= Y] f d^Xi /d 3 X 2 ...d 3 X m ^ Al ,a 1 ,...,A m ,a m (xi,...,X m ) 

Ai ,ai ,. . .,A m ,a m 

X(^l^ 1 ,a 1 (xi)..-^ m , am (x m ) |0) 

/ j rf3x 2---^a;m.|^A 1 ,a 1 ,...,A„ i ,a m (Xl,...,X m )| 2 

Ai ,ai ,. . .,A m ,a m ^ 

wl. (115) 

The approximation in the last line stems from the fact that most of the support of 
the wavefunction <£>Ai,ai,...,A m ,a m (xi, . . . ,x m ) is within the region B m C M 3m . For similar 
reasons we have 

<0|C„(B)&(B)|0)«1. (116) 
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We also have 



Ai,ai,...,A m 

(m-1) 



x 



m 



"V- , Ai,ai,...,A m ,a m ( x l; • ■ • j x m) 



+ — / rf 3 x5 {A) (x-Xi)^ Aliaii ... iAmiam (x,x 2 ...,x m ) ] 

m J B J 



m 



m 



m 



(117) 



Here the approximation stems not only from the fact that most of the support of the 
wavefunction is within B m , but also from the fact that most of the support of the 
function 5^ A ' is located within a region of the order of the Planck length 1/A, so that 
we can approximate the expression by replacing it with the ordinary (^-function. 
We further have that 

&\F d (B)Cl(B)\0) = (0|3„F d (£)qt(£)|0) 

= (0| [d v ,F d {B)] Cl(B)\0) + (0\F d (B)C v Cl(B)\0) 

= m(0\C^B)Cl(B)\0) + ( V \Cl(B)\0)(0\F d (B)\0) 



m + uq(B) 



'im 



where in the last line we respectively used (116) and (115). With the help of this relation 
and (117) we have that 

(y\CUB)F d {B)\Q) = (cp\F d (B)Ci(B)\0) - (<p\ \F d {B),&JB)] |0> 



= (<p\F d {B)6l{B)\0)-m(<p\C ip {B)\0) 
m + uq{B) — m = no(B) . 

C.3 Expectation value and standard deviation 

Using the approximations, we obtain the expectation value 

(<p\F d {B)\<p) = ( V \F d (B)dl\0) 

= (<p\ [F d (B),Cl] \0) + (<p\d},F d (B) |0) 
= m(v\Cl(B)\0)+n (B) 

~ Tfl -\- Uq{B) 



;ii9) 



(120) 
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and with 



(< P \F d (B) 2 \^} = ( l p\F d (B) 2 dl\0} 

= (<p\F d (B) [F d (B),dl] |0> + ( V \F d (B)dlF d (B) |0> 

= m(tp\F d (B)Cl(B)\0) + (<p\ [F d (B), 5+] F d (B) |0> + (<p\ C^Bf |0> 

= m{ V \F d {B)& v {B)\Q) + m (<p\ Cl(B)F d (B)\0) + (0| F d (B) 2 |0> 

« m 2 + 2mn (B) + (0\F d (B) 2 |0) , (121) 

we find the standard deviation 



A v F d (B) = (<p\ F d {Bf \tp) - (ip\ F d (B) \<p} 2 « A F d (B) . (122) 

Completely analogously one can calculate the expectation value (ip\ F d (B) \<p) and 
standard deviation A^F^B) for a region B C M 3 \ B. The result is: 

(if\ F d (B) \<p) « ^(5) , A^(fl) « AoF d (5) . (123) 
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